Analysis of matter suppression in collective neutrino oscillations 
during the supernova accretion phase 
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The usual description of self-induced neutrino flavor conversions in core collapse supernovae (SNe) 
is based on the dominance of the neutrino density over the net electron density n^. However, this 
condition is not met during the post-bounce accretion phase, when the dense matter in a SN is piled 
up above the neutrinosphere. As recently pointed-out, a dominant matter term in the anisotropic 
SN environment would dephase the flavor evolution for neutrinos traveling on different trajectories, 
challenging the occurrence of the collective behavior in the dense neutrino gas. Using the results 
from recent long term simulations of core-collapse SN explosions, based on three flavor Boltzmann 
neutrino transport in spherical symmetry, we find that both the situations of complete matter 
suppression (when n^) and matter- induced decoherence (when ?ie f^njj) of flavor conversions 

are realized during the accretion phase. The matter suppression at high densities prevents any 
possible impact of the neutrino oscillations on the neutrino heating and hence on the dynamics 
of the explosion. Furthermore, it changes the interpretation of the Earth matter effect on the SN 
neutrino signal during the accretion phase, allowing the possibility of the neutrino mass hierarchy 
discrimination at not too small values of the leptonic mixing angle On (i.e. sin^ 6*13 > 10^^) . 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

Massive star explosions are an active subject of re- 
search, in terms of core-collapse supernova (SN) models. 
The supernova problem is related to the revival of the 
stalled bounce shock, which forms when the collapsing 
stellar core reaches nuclear matter density and bounces 
back. It propagates out of the core and stalls on its way 
out due to the continuous energy loses via neutrino emis- 
sion and dissociation of heavy nuclei. Independently from 
the explosion mechanism (see, e.g. [H^) the total energy 
emitted in neutrinos (v) and antineutrinos [v) during a 
SN is of the order of several 10^^ erg, which makes a SN 
the most powerful neutrino source in the Universe. The 
total duration of the neutrino emission, known as the 
neutrino burst, can last up to 30 seconds. This v signal 
represents a powerful tool to probe fundamental neutrino 
properties as well as the dynamics of the explosion . 

The role of astrophysical messengers played by neutri- 
nos during a stellar collapse is largely associated with the 
signatures imprinted on the observable SN neutrino burst 
by weak processes and flavor conversions occurring deep 
inside the star. In the last few years, it has been under- 
stood that the description of neutrino flavor conversions 
in supernovae, based on the only Mikheyev-Smirnov- 
Wolfenstein (MSW) effect with the ordinary matter 
was incomplete. In particular, in the deepest supernova 
regions the neutrino density is so high that the neutrino- 
neutrino interactions dominate the flavor changes. Even 
if the important role played by the neutrino-neutrino 
interactions for the flavor evolution in supernovae was 
pointed out since a long time ago [lol - [l^ . only recently 
the first numerical experiments have been performed to 



characterize realistically these effects [l^ [T^ . These sem- 
inal investigations have stimulated a still streaming tor- 
rent of activities (see, e.g. [13, for a recent review). The 
general result of these studies is that neutrino-neutrino 
interactions create a large potential for the neutrinos, 
which causes large and rapid conversions between differ- 
ent flavors. The transitions occur collectively, i.e. in a 
coherent fashion, over the entire energy range. The most 
important observational consequence of v-v interactions 
is a swap of the Vf, and Ve spectra with the non-electron 
Vx and Vx spectra in certain energy ranges (l6l - [3^ . 

The basic idea behind the self-induced neutrino oscilla- 
tions in supernovae is that these can produce significant 
flavor conversions close to the v decoupling region (i.e., 
at small radii) even in the presence of a large matter 
density [Sq . This would have inhibited any possible 
flavor transition, far from the MSW resonances (i.e., far 
outside the neutrino-matter decoupling region) at larger 
radii (36| . Subsequently, it has been realized in 37 1, and 
then conflrmed in [s^ , that this original idea was only in 
part true since the matter density cannot be arbitrarily 
large before it affects collective flavor conversions after 
all. Indeed, for the strongly anisotropic emission of neu- 
trinos streaming-off the SN core, the matter effect is tra- 
jectory dependent, inducing different oscillation phases 
for neutrinos traveling in different directions. When 
the neutrino-neutrino interaction is sufficiently strong, it 
forces neutrinos reaching the same position from different 
paths to have the same oscillation phase, overcoming the 
trajectory-dependent dispersion induced by the matter. 
Conversely, when the electron density Tig significantly ex- 
ceeds the neutrino density rii, , it has been shown that the 
large phase dispersion induced by the matter would sup- 
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press the collective phenomena. Depending on the elec- 
tron density, the matter suppression can be complete, 
when Tie S> ?T.jy, or partial when the matter domination is 
less pronounced. Finally, when ne>7iy, the interference 
of the two comparable effects would lead to the decoher- 
ence of the collective neutrino flavor changes, producing 
an equal mixture between the oscillating electron and 
non-electron neutrino species (STj . 

Most of the numerical simulations of the self-induced 
flavor conversions in SN explosions have been focused 
on the late times cooling phase (at post-bounce times 
tpt > 1 s) , when the matter potential is smaller than the 
neutrino-neutrino potential. In such a situation, collec- 
tive oscillations would develop without any matter hin- 
drance [l^, [l^, [311, • During the earlier phase before 
the onset of an explosion, which is determined by mass 
accretion (tph<0.5 s), the electron density is not negli- 
gible with respect to the neutrino density, as one can in- 
fer from the following back-of-the-envelope calculation. 
Given the mass accretion rate 



dm 2 
m = -^v = Airpr v 
or 



(1) 



in terms of the the rest-mass density p at radius r, and 
of the matter velocity v, one gets 



cm3 VO-4 Mq/s 



10^ km^ /lO^ km/s 



(2) 

for typical values in the region ahead the shock-wave. 
Then, the net electron density at r ~ 10^ km can be 
estimated as 



pYe/niB 10"^ cm 



(3) 



where the electron fraction Ye ~ 0.5 and tub is the nu- 
cleon mass. On the other hand, during the accretion 
phase, where the charged-current reactions dominate, the 
far-distance electron (anti)ncutrino luminosities can be 
expressed via the change of the gravitational potential at 
the neutrinosphere ^ with radius Tj, as follows. 
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at r ~ 10^ km. From this simple estimation, can be 
smaller or on the same order than and hence we cannot 
a-priori assume that rie <C holds during the post- 
bounce accretion phase. This result is confirmed by many 
different SN simulations [39l - l42j . It is a generic feature 
that applies generally to core collapse SNe of massive 
iron-core progenitors. As an important consequence, it 
is not guaranteed a priori that the results of the neutrino 
flavor evolution obtained during the cooling phase, could 
be directly applied to the accretion. 

The post-bounce accretion phase, where the neutrino 
fluxes are large, represents a particularly interesting sce- 
nario for detecting signatures of neutrino flavor trans- 
formations. Indeed, also the flavor-dependent flux dif- 
ferences are large with a robust hierarchy for the neu- 
trino number fluxes, F^^ > F-^^ 3> i^i^^, where Vx in- 
dicates the non-electron flavors. In such a situation, if 
dense matter effects are neglected and the flavor asym- 
metries are enough to prevent multi-angle effects [43| . 
the self-induced flavor oscillations seemed to have a clear 
outcome, producing a complete exchange of the elec- 
tron and non-electron flavors for almost all antineutrinos 
in the inverted neutrino mass hierarchy (IH, Am^jj,-, = 
rn^ — m\ 2 < 0) £^nd a s pec tral split in the energy distribu- 



tions of the neutrinos jl8| . In the normal neutrino mass 
hierarchy (NH, Am^^j^^ > 0), they would leave the neu- 
trino spectra unaffected. How this picture would change 
if the effect of the dense matter is included, is subject 
of investigation of the present article. It develops and 
clarifies the results recently presented in our [3] . 

We take as benchmark the results of the recent long- 
term SN simulations from Ref. [3l , to characterize the SN 
neutrino signal as well as the matter density profile evolu- 
tion. We consider three cases corresponding to different 
supernova progenitor masses, the 8.8 Mq 0-Ne-Mg-core 
and the 10.8 M© and 18.0 M© iron-cores. For these three 
cases wc find that the matter density is always larger than 
or comparable to the neutrino density, during the accre- 
tion phase. This finding suggests the possible hindrance 
(or decoherence) of collective neutrino flavor transitions 
for matter densities higher than (or comparable with) the 
neutrino densities during the accretion phase. 



10^ 



s Mq) VO.4 Mq/s /V ^ 



10^ km 



(4) 



Therefore, the neutrino particle flux F^^^p^ = 
Lp^,i,J{Ep^,uJ ~ 10^^ s-\ for {Ep^,p^) = 11 MeV, im- 
plying 



47^7-2 



lO'^'cm 



(5) 



^ The neutrinospheres are the neutrino energy and flavor depen- 
dent spheres of last scattering. 



The plan of our paper is the following. In Section II we 
describe the original supernova neutrino emission during 
the accretion phase and in Section HI we introduce the 
schematic supernova neutrino model we use to charac- 
terize the v emission geometry, the neutrino potentials 
and the different radial ranges where self-induced oscil- 
lations and matter effects arc relevant. The results of 
our analysis on matter suppression for the three different 
SN progenitor masses are presented in Section IV and 
the impact of the dense matter suppression on the inter- 
pretation of the Earth matter effect on the SN neutrino 
burst is discussed in Section V. Finally, we close with 
comments of our results and the conclude in Section VI. 



3 



\ 
cn 

0) 

o 



7^ 



8.8 Msun 




^ 16 
^ 14 
A 12 

UJ 10 

V „ 



0,05 0.1 0.15 0,2 0,25 



0,05 0,1 0,15 0.2 0,25 




0,1 0,2 0,3 0,4 0.5 0,6 
tpb (s) 



0,1 0,2 0,3 0,4 0,5 0,6 
tpb (s) 



FIG. 1: Neutrino luminosities (left panels) and average en- 
ergies (right panels) for the core-collapse SN simulations con- 
sidered, where we distinguish (continuous line), Ve (dashed 
line) and (dotted line). The three cases are shown corre- 
sponding to different values of the supernova progenitor mass: 
8.8 Mq (top), 10.8 Mq (middle) and 18.0 M© (bottom). 



II. NEUTRINO SIGNAL FROM 
CORE-COLLAPSE SUPERNOVA EXPLOSIONS 

We investigate core-collapse supernova simulations of 
massive progenitor stars of 8.8, 10.8 and 18 M©. The 
first one belongs to the class of 0-Ne-Mg-core progeni- 
tors [i^ and represents the threshold between ther- 
monuclear explosions and core-collapse supernovae [1, 0] . 
The latter two are iron-core progenitors |47[ . All models 
were evolved consistently through core collapse, bounce 
and the early post-bounce phase up to several seconds 
after the onset of explosion [3]. The core-collapse model 
is based on general relativistic radiation hydrodynam- 
ics that employs three flavor Boltzmann neutrino trans- 
port in spherical symmetry and a sophisticated equation 
of state for hot and dense nuclear matte r l48t (for de- 
tails about the supernova model, see 0, |49{ and refer- 
ences therein). Explosions of massive iron-core progen- 
itors cannot be obtained in spherically symmetric su- 
pernova models. In order to trigger the explosions for 
the 10.8 and 18 Mq progenitor models, the heating rates 
are artificially enhanced in the gain region where neutri- 
nos deposit energy in order to revive the stalled bounce 
shock. 

Fig. [T] shows the evolution of the neutrino luminosities 



(left panels) and average energies {E^^ (right pan- 
els) , during the post-bounce accretion phase for all mod- 
els under investigation. Here, = (i^e, t^e, '^k) where Vx 
indicates both (/i, T)-neutrinos and antineutrinos. The 
large electron fiavor neutrino luminosities 0(10^^) erg/s 
are due to the continuous mass accretion at the neutri- 
nosphercs, where the electron fiavor neutrino luminosities 
are dominated by the charged current reactions. The 
slightly larger luminosity over the v^^ luminosity and 
the magnitude of the differences obtained in the mod- 
els discussed here, is an active subject of research and 
may change slightly applying improved weak rates and 
multi-dimensional supernova models. The increasing and 
decreasing electron flavor luminosities reflect the prop- 
agating bounce shock (and the mass accretion rate at 
the neutrinosphercs) , which contracts and expands ac- 
cordingly during the accretion phase driven by neutrino 
heating/cooling. Since muonic charged current processes 
are suppressed, i'^/t-'s are produced only after bounce 
via pair-processes and the (/i/T)-neutrino luminosities 
are generally smaller than the electron flavor neutrino 
luminosities. For the average energies, the following hier- 
archy holds (Ev^) < {Ejy^) < {EyJ) for all models under 
investigation during the accretion phase. Furthermore, 
the average energies rise continuously for all neutrino fla- 
vors during the accretion phase illustrated in Fig. [1] (right 
panels). 



After the onset of explosion, when mass accretion van- 
ishes, the situation changes. The luminosities and aver- 
age energies of all flavors decrease on a longer timescale 
on the order of seconds. For the 0-Ne-Mg-core, we 
consider the neutrino signal only up to tpb=0.25 s (top 
panel of Fig. [T]). Note that mass accretion at the neu- 
trinosphercs vanishes already at about fpb=0.03 s. This 
determines the onset of the explosion for this low-mass 
progenitor. For the 10.8 and 18 Mq progenitor mod- 
els discussed here, we show the neutrino signal up to 
ipb=0.6 s. For these cases, the onsets of explosion occur 
at about <pb=0.36 s, due to the more massive envelopes 
surrounding the iron-core that lead to a more extended 
accretion (and hence neutrino heating) phase. Note the 
sharp drop of the luminosities and average energies of 
all neutrino flavors in Fig. [1] after the onsets of explo- 
sion. These are due to the sudden flip of matter velocities 
from infall to expansion when the explosion shock passes 
through the distance of 500 km, where the observables 
are measured in a co-moving frame of reference. We note 
further that these results suggest lower average energies 
than often assumed in the literature [50| and a less pro- 
nounced spectral hierarchy, in particular during the later 
proto-neutron star cooling phase. The results obtained 
for the low mass 0-Ne-Mg-core collapse supernova ex- 
plosion and long term simulation of the proto-neutron 
star cooling phase, arc in qualitative and quantitative 
good agreement with recent simulations performed by 
the Garching-group [sij . 
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FIG. 2: 10.8 Mq progenitor mass. Angular distributions 
for Ve. at tpb=0.4 s for Srms = 11-5 MeV at different radii: 
r=20 km (dotted line), r=50 km (continuous line), r=70 km 
(dashed line). 



III. SCHEMATIC SUPERNOVA MODEL 

We present here our schematic assumptions to charac- 
terize the flavor conversions during the accretion phase. 
In particular, we describe our model for the SN neutrino 
emission geometry, the neutrino number fluxes of differ- 
ent flavors, the neutrino- neutrino and matter potentials 
and the different oscillation regimes. 

A. Neutrino emission geometry 

We assume a spherically symmetrical source ("SN 
core") that emits neutrinos and antineutrinos like a 
blackbody surface, from a neutrinosphere at radius r = 
Tjy. We define to be the radius where the neutrino 
radiation field is assumed to be "half-isotropic" , i.e. all 
outward- moving angular modes are equally occupied [43| . 
We remind that this conventional definition of the neu- 
trinosphere, used in the context of the neutrino flavor 
conversions, is only intended to fix a boundary condi- 
tion for the subsequent fiavor evolution. We assume that 
weak processes take place at higher densities inside the 
neutrino-matter decoupling sphere and are negligible out- 
side the decoupling sphere where neutrino oscillations can 
occur. The half-isotropic definition for the decoupling 
spheres does not necessarily coincide with the definition 
of the neutrinosphere as neutrino last scattering surface, 
defined as the radius where the optical depth becomes 
2/3. However, in the context of collective oscillations 
which start at radii much larger than the assumed bound- 
ary, neutrinos are safely in a free-streaming regime (3l| . 
Therefore, we do not have to worry about the details of 
the V decoupling. 

Instead of choosing arbitrarily the neutrinosphere ra- 



dius, we fix it consistently with the SN simulations which 
provide the angular distributions of the different neu- 
trino species as a function of time and energy at differ- 
ent radii. For definiteness, we have taken the z^e's dis- 
tribution as representative for all the different flavors, 
since electron neutrinos reach the free-streaming regime 
at larger radii compared to the other neutrino species 
during the accretion phase. The angular distributions 
are functions of the energy. For definiteness, we con- 
sider as reference the root mean square (rms) energy rel- 
evant for the average neutrino-nucleon interaction rates, 
that determine the neutrino angular distributions. As an 
example, in Fig. [2] we plot the angular distribution 
f {cos 9r) — dn^ /d cos 6r, where dr is the zenith angle 
of a given mode relative to the radial direction at dis- 
tance r. We consider E'lms = 11-5 MeV at three different 
radii at the post-bounce time tpb — 0.4 s, taken form the 
10.8 M0 model. As expected, at small radii (r — 20 km 
in the Fig. [5]) the angular distribution is isotropic since 
neutrinos are in a trapping regime. They are isotropically 
emitted in all directions. At large radii (r = 70 km in the 
Fig. HI), neutrinos are free-streaming and their angular 
distribution becomes forward peaked. We schematically 
assume as neutrinosphere radius, the one at which the 
I'e's angular distribution has no longer significant back- 
ward flux, i.e. a few % of the total one (at r = 50 km 
in the Fig. [5]). Even if there the real angular distribu- 
tion is not half-isotropic, we assume it to characterize 
the further fiavor evolution. For our purpose of evaluat- 
ing the impact of the matter effects on the self-induced 
fiavor conversions, this simplified choice is conservative. 
Indeed, the real angular distributions, which are more 
forward-peaked than what we are assuming, would re- 
duce the real strength of the neutrino-neutrino potential 
with respect to what we are assuming. 

We stress that this procedure for fixing the neutri- 
nosphere has to be taken only as an empirical prescrip- 
tion. However, we performed several neutrino flavor os- 
cillation analysis changing by a factor of three the neu- 
trinosphere radius. We find that the results on matter 
effects remained basically unchanged, indicating a weak 
dependence of the oscillation analysis from the decou- 
pling region. 



B. Neutrino number fluxes for different flavors 

The neutrino particle flux for the different flavors is 
defined as follows F^^ = L,^^/{E^^). In order to charac- 
terize the flux asymmetries among the different v species 
we introduce the parameter [i^l 



F^^ 



Fv 



(6) 



where we assume F^^ — F-^^ . 

In Fig. [3] we show the neutrino number flux F^^ (left 
panels) and asymmetry parameter e (right panels) for the 
three benchmark SN simulations based on the different 
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FIG. 3: Neutrino particle fluxes (left panels) and flavor asym- 
metry parameter e (right panels) for v^. (continuous line), 
(dashed line), (dotted line). 



progenitor models. The particle fluxes present the ex- 
pected hierarchy, Fy^ > F^^ :s> Fi,^ , during the accretion 
phase. The first part of the hierarchy is caused by the 
deleptonization of the collapsed stellar core whereas the 
second is due to the absence of charged-current interac- 
tions for neutrino species other than v^. and V,,- As a con- 
sequence of this hierarchy, the electron and non-electron 
neutrino particle fluxes differ by almost a factor of two 
during the accretion phase. Conversely, flux differences 
tend to become much smaller during the cooling of the 
SN remnant. 

Passing to the asymmetry parameter e [Eq. (jS])], the 
strong excess of Ve 's during the early post-bounce delep- 
tonization (<pb<0.05 s), increases e 1. Then, e drops 
and reaches values between 0.3-0.5 during the accretion 
phase (see Fig. [3]). After that, since the flux difference 
F^^ — F-^^ drops more rapidly than F^^ — F^^ , e rises again 
and becomes larger than 1 for the iron-core SNe. We will 
see that this behavior will have important consequences 
on the development of the self-induced flavor conversions. 



C. Non-linear neutrino flavor mixing 



flavor effects, associated with the solar sector, are small 
for the neutrino flux ordering expected during the accre- 
tion phase [H, nil . 

In the normal mass hierarchy and for the spectral or- 
dering of the accretion phase, no self-induced flavor con- 
version will occur and thus dense matter cannot produce 
any sizable new effect. Therefore, in the following we will 
always refer to the inverted mass hierarchy case. 

The impact of the non-isotropic nature of neutrino 
emission on the self-induced flavor conversions is taken 
into account by "multi- angle" simulations jJLfij, where one 
follows a large number [0(10^)] of neutrino trajectories. 
The ^'s emitted from a SN core naturally have a broad 
energy distribution. However, this is largely irrelevant 
for our purposes, since large matter effects would lock 
togheter the different neutrino energy modes, both in 
case of supression [s^l and of decoherence [i^ of collec- 
tive oscillations. Therefore, to simplify the complexity of 
the numerical simulations, we assume all i^'s to be repre- 
sented by a single energy, that we take E — lb MeV. This 
approximation is reasonable since our main purpose is to 
determine if the dense matter affects the development of 
the self-induced transformations. This choice results in 
the neutrino vacuum oscillation frequency 



Ato„ 



2E 



= 0.4 km 



(7) 



The strength of the neutrino-neutrino interaction, nor- 
malized at the neutrinosphere, is parametrized by |43| 



Hr = \/2Gf [nv^ {r) - n-j^^ (r)] 



= 7 X 10^ km"^ 



15 MeV / 10 km 



lO^^erg 



(8) 



where n^^ (r) is the flux of the neutrino species Va at 
radius r (for a definition of Uj,, see Eq.(l5]))- The matter 
potential is represented by [36[ 



A,. = V2GFne{r) = 1.9 x 10^ km^^ 



0.5 J \10^° g/cm3 



(9) 



encoding the net = n^- — 71^,+ electron density, where 
Ye — Yg- — y"g+ is the net electron fraction and p is the 
matter density. The neutrino term /x^ declines always as 
due to the geometric dilution outside the decoupling 
spheres, whereas is given by the detailed matter profile 
from the SN simulations. 



Our description of the non-linear neutrino flavor 
conversions is based on the two-flavor oscillation sce- 
nario, driven by the atmospheric mass-square difference 
Am^tjjj — 2.6 X 0"'^ eV^ and by a small (matter sup- 
pressed) in- medium mixing OcS = 10~^ [s^ [ssj. Three- 



D. Oscillation regimes 

In the absence of matter suppression, collective neu- 
trino flavor transformations will start outside the syn- 
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FIG. 4: 10.8 M© progenitor mass. Time evolution of (con- 
tinuous line), Tsync (dashed curve) and Tend (dotted curve). 



chronization radius, given by [43 



1/2 



1/4 



(10) 



and are expected to develop at least until radii when the 
neutrino-neutrino interaction strength becomes compa- 
rable to the vacuum term, i.e. at llSlI 



1/4 



^end 

r,. 



2uj 



(11) 



We remind that outside rgync the multi-angle nature 
of the neutrino trajectories in a supernova can lead to 
self-induced flavor decoherence between different angular 
modes. However, it has been shown in [4^ that if the 
asymmetry parameter e is significantly large, multi-angle 
effects are suppressed and self-induced neutrino oscilla- 
tions exhibit a collective behavior. We will come back to 
this point in more details when commenting our numer- 
ical results. 

Using the above definitions, in Fig. |4] we present 
the neutrinosphere radius , the synchronization radius 
?'sync and the radius Tend at which collective effects sat- 
urate, for the case of the 10.8 SN simulation. We 
see that ^ 10^ km during the accretion phase and 
drops to ~ 30 km at the beginning of the cooling phase. 
This contraction of the neutrinosphere radius reflects the 
gradual shift of the v opacities. These change from be- 
ing absorption dominated during the accretion, to being 
scattering dominated during the cooling [53, [5a| ■ 

Collective oscillations are expected in the range r E 
[?'sync, '"end]- This range is at r ~ [600,1500] km at 
^pb = 0.1 s when the neutrinosphere radius > 10^ km, 
pushing to larger radii the flavor conversions. Then, fol- 
lowing the contraction of the neutrinosphere radius, the 
range of conversions shifts at smaller radii, moving to 



r - [200, 500] km at tp^ ~ 0.4 s. Finally at tpb ^ 0.6 s, 
the lower neutrino luminosity and the larger asymmetry 
parameter e both conspire to push rgync towards rend [see 



Eqs. (Unil-dlll)]- 

In the range 



' sync ; 



'"cndj 



the self-induced neutrino fla- 



vor conversions can be affected by the dense matter 
when llll 



rip > n-r 



(12) 



In particular, when the two densities are comparable 
matter effects induce multi-angle decoherence among dif- 
ferent neutrino modes, leading to a flavor equilibration 
among the different species. When the net electron den- 
sity is significantly larger than the neutrino density, col- 
lective oscillations are suppressed at all [37| . 

We remark that previous numerical studies of self- 
induced flavor conversions typically assumed smaller neu- 
trinosphere radii, namely = 0(10 km) (see, e.g., [43|). 
This choice is more appropriate for the cooling rather 
than for the accretion phase. As a consequence of this dif- 
ferent input, in the current work flavor conversions start 
at larger radii than typically assumed in previous studies. 



IV. MATTER EFFECTS FOR DIFFERENT 
PROGENITOR MASSES 

In this Section we will present our results concerning 
the effects of dense matter on of the self-induced neutrino 
oscillations for the three representative SN simulations 
based on the different progenitor masses. 



A. 10.8 Mq progenitor mass 

We start our investigation with the case of the 10.8 M0 
iron-core SN. In Fig. [S] we show the net electron den- 
sity rie (left panel) and the difference of neutrino den- 
sities n-p^ — n-jy^ (right panel) entering in the potential 
fir [Eq- for different post-bounce times. While the 
decline of the neutrino density is r~^, the electron den- 
sity presents a more complicated behavior. In particular, 
one can recognize the abrupt discontinuity in rig asso- 
ciated with the supernova shock-front that propagates 
in time. The shock-wave initially dissipates its energy 
due to heavy nuclei dissociation and eventually stalls at 
tph — 0.05 s. The standing accretion shock is revived via 
neutrino heating, on a timescale on the order of 0.1 sec- 
onds, and expands accordingly as well as contracts via 
neutrino cooling. After 0.35 seconds post-bounce, neu- 
trino heating succeeds and the standing accretion shock 
turns into a dynamic shock with positive matter veloci- 
ties. It initiates the onset of explosion at about 0.4 sec- 
onds post-bounce. 

For this SN model we find that during the accre- 
tion phase the matter density in the post-shock region 
declines slower than the neutrino density, typically as 
~ j.-i-5^ From the comparison between the electron and 
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FIG. 6: 10.8 M0 progenitor mass. Radial evolution of the ratio R between electron and neutrino densities at different 
post-bounce times. The two dashed vertical strips delimit the position of rsync (left line) and Tend (right line). 
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the neutrino densities, we realize that at the different 
post-bounce times considered, is always larger than 
or comparable to nv^ — riv-^ ■ It suggests that one cannot 
ignore matter effects on self-induced flavor transforma- 
tions during the accretion phase. 

In order to quantify the relative strength of the electron 
and neutrino densities, in Fig. [S] we show the ratio 



R 



n-v. - m 



(13) 



as a function of the radial coordinate r at different 
post-bounce times for ipb G [0.1,0.6] s. The range 
''sync < r < Tend IS delimited with two vertical dashed 
lines. The values of rgync and rend determine the possi- 
ble range for the self-induced flavor conversions and the 
shock radius Tgh denotes the abrupt drop in the electron 
density. Therefore, their relative position is crucial to as- 
sess the impact of matter effects. In the expected oscilla- 
tion range, i? 3> 1 will imply a strong matter dominance 
in the flavor conversions and thus complete suppression 
of the self-induced effects. Instead, when electron and 
neutrino densities are comparable (i?>l), decoherence 
will occur for the collective oscillations. 

The ratio i?, being very large behind the shock front, 
prevents flavor conversions in this region. However, the 
ratio can go down to i? > 1 for r > Tgh, leading to matter- 
induced decoherence and thus partial flavor changes. 

Let us discuss in more detail what occurs at differ- 
ent post-bounce time snapshots in Fig. [51 At very early 
times (tpb — 0.1 s) the matter term is strongly domi- 
nant also behind the shock- front (i? 3> 1). Under these 
conditions oscillations are always blocked. Then, at in- 
termediate times (tpb = 0.225,0.3 s) the matter den- 
sity in the post-shock region, where flavor conversions 
are possible, is dropping faster than the neutrino one. 
Therefore, the ratio R drops at 1-2 in this range and 
matter-induced decoherence is possible in this case. Sub- 
sequently, at ipb = 0.325 s, oscillations are suppressed 
behind the shock front, but then decoherence will de- 
velop at larger radii (r>300 km) when -R>1. Eventu- 
ally at later times (tpb = 0.4,0.6 s), since the shock has 
resumed its forward motion, the region relevant for the 
oscillations is at r < Tgh, where i? 3> 1. In this situ- 
ation self-induced oscillations will be suppressed. From 
these different snapshots we realize that R has a peculiar 
non-monotonic behavior as a function of time. It sug- 
gests a time-dependent pattern for the matter effects on 
the self-induced transitions during the accretion phase, 
namely complete-partial-complete suppression. 

In order to conflrm these expectations, we have per- 
formed a multi-angle numerical study of the equations 
for the neutrino flavor evolution in the schematic model 
described in Section III. Our treatment closely follows 
the one presented in Ref. [S^l to which we address the 
interested reader for further details. We only mention 
here that in order to achieve convergence in our simu- 
lations we had to simulate lO'^ neutrino angular modes. 
In Fig. [7] we show the radial evolution of the X/g survival 
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FIG. 7: 10.8 Mq progenitor mass. Radial evolution of the 
survival probability Pee for electron antineutrinos at diflferent 
post-bounce times for the multi-angle evolution in presence 
of matter eflFects (continuous curve) and for He = (dashed 
curve). 



probability Pee for different post-bounce times, obtained 
taking into account the effects of the SN matter profile 
(continuous curve). For comparison, we show also the 
results obtained setting Ue = (dashed curve). 

In the case with Ue = 0, for the given flavor asym- 
metry e > 0.3 we would have expected the "quasi-single 
angle" behavior described in Ref. [i^, where after the 
onset of the conversions at r — rgync , the survival prob- 
ability Pee declines smoothly approaching zero at large 
radii. However, in the situation we are studying flavor 
conversions develop at radii larger than what is typi- 
cally shown in previous works (see, e.g., 43!]). There- 
fore, the evolution is more adiabatic (i.e. the evolu- 



tion length scale 



[3d]). As a consequence, ef- 



fects of self-induced multi-angle decoherence have more 
chances to develop in this case, producing some small 
disturbance in the smooth decline of the survival proba- 
bility at large radii (visible at tpb — 0.1,0.3,0.325 s for 
r > 700 km) . This flnding is potentially interesting, how- 
ever, since matter effects will anyhow dramatically alter 
this picture we have not performed a systematic study 
on this (sub-leading) self-induced decoherence. 

Passing now to the matter case, we see that at 
tph = 0.1,0.4,0.6 s the flavor oscillations are completely 
blocked, since i? ^ 1 in the conversion range. For the 
other three intermediate times (tpb = 0.225, 0.3, 0.325 s), 
the presence of a large matter term at rgync will signifl- 
cantly delay the onset of the flavor conversions with re- 
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tpb(s) 

FIG. 8: 10.8 Mq progenitor mass. Upper panel: Time evo- 
lution of the ratio i?sync (dashed (continuous 
curve) and i?cnd (dotted curve). Lower panel: Time evolution 
of the radial position rsync (dashed curve) , rmin (continuous 
curve), Tend (dotted curve) and r-shock (dot-dashed curve). 



sition of the minimum value i?min is at 



^sh if ^sync ^ ^sh ^ ^cnd 



' sync 



otherwise 



(15) 



For tpb < 0.3 s, since the neutrino density grows and the 
electron density decreases in time, i?min decreases mono- 
tonically from ~20 at tph = 0.05 s to ~l-2 at tph — 0.3 s 
post-bounce, changing from complete to partial matter 
suppression. Then, since at the end of the accretion 
phase the neutrino density decreases faster than the elec- 
tron density, i?min rises again to ~20 at tph — 0.4- 
0.6 s, suppressing flavor conversions. The behavior of 
Rmin shows in a compact way the succession of phases 
with complete-partial-complete matter suppression in 
the time evolution of the neutrino self-induced conver- 
sions during the accretion phase. 

A similar analysis to the one presented above based on 
Fig. [51 may allow to obtain a first answer about the role 
of the matter on the self-induced effects. This schematic 
approach does not require a detailed numerical study of 
the neutrino flavor evolution. For this reason, we find 
it useful to compare the role of matter in different SN 
simulations based on different progenitor models. In the 
following, we will apply such analysis to the other SN 
simulations from Ref. J,] based on different progenitor 
masses. 



B. 18.0 Mq progenitor mass 



spect to the case with rig = 0. Then, at larger radii 
(r > 700 km) when R > 1-2, matter effects produce a par- 
tial suppression of the flavor conversions, with a tendency 
toward flavor decoherence. When the latter is complete, 
it leads to Pee — >■ 1/2 (at tpb — 0.3,0.325 s), implying a 
complete mixture between Ve and . 

Finally, we would summarize our results with a figure 
of merit that captures in a compact way the interplay 
between matter and self-induced effects as a function of 
time. Therefore, in Fig. [5] we show the ratio R (upper 
panel) at its minimum value 

Rmin = min , (14) 

in the range r S [j'sync; ''end] relevant for collective con- 
versions. This value provides a conservative estimation 
for the impact of matter effects. Indeed, if i?min ^ 1 one 
should expect a complete matter suppression of the flavor 
changes. Conversely, when Rmin ^ 1 the matter suppres- 
sion will be partial and will lead to flavor decoherence. 
For comparison we also show the values of R at the two 
radii rgync and rend- Moreover, in the lower panel we 
represent the time evolution of the radial positions rgync, 

''end) '^min and Tgh- 

One realizes that R presents a non-trivial behavior in 
the range between Vgync and rend- In particular, the po- 



As further example, we consider the SN simulation of 
the more massive 18 Mq iron-core progenitor. In Fig. [5] 
we show the net electron density rig (left panel) and the 
difference of neutrino densities n-j^^ — n-^^ for different 
post-bounce times. From the comparison with Fig. [5] 
we see that the time evolution of the electron and neu- 
trino density profiles is similar to the case 10.8 Mq SN 
simulation. Therefore, we expect a similar impact for 
the matter effects on the neutrino flavor evolution. In 
Fig. \W\ we show the time evolution of i?min in the same 
format as in Fig. |8l We see that the behavior of i?min 
is also very similar to the previous case described above. 
Therefore, one expect an analogous pattern of complete 
0.05 <tpb< 0.2 s), partial (0.2 < ipb < 0.35 s) and com- 
plete (0.35 <tpb < 0.6 s) matter suppression in the self- 
induced flavor conversions. 

We have repeated the same analysis also for a 15.0 Mq 
iron-core progenitor mass, giving a not-exploding super- 
nova. In the time range before the recoUapse of the star 
(ipb<0.3 s) the results obtained are similar to the one 
shown for the two exploding cases. Therefore, for the 
sake of the brevity, we will not show them here. 

C. 8.8 Mq progenitor mass 

Finally, we pass to analyze the case of the SN of the 
low mass 8.8 Mq 0-Ne-Mg core. In Fig. [TT]we present 



10 



I 



o 



CD 
C 



-1 
1 






-2 
10 




IN. ~E 


-4 
10 






-6 
10 






10^ 


r (km) 




r (km) 



FIG. 9: 18.0 M0 progenitor mass. Radial evolution of the net electron density Ue (left panel) and of the neutrino density 
difference n-^ — nj;^ (right panel) at different post-bounce times. 
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FIG. 10: 18.0 Mq progenitor mass. Upper panel: Time 
evolution of the ratio Rsync (dashed curve), i?min (continuous 
curve) and i?cnd (dotted curve). Lower panel: Time evolution 
of the radial position rgync (dashed curve) , Tmin (continuous 
curve), Tend (dotted curve) and rshock (dash-dotted curve). 



the net electron density rig (left panel) and the difference 
of neutrino densities n-jy^ — n-jj^ (right panel) for different 
post-bounce times. We realize that since in this case 
the matter density of the envelope is very low compared 
to the iron-core progenitors, the electron density profile 
above the core is very steep, declining as ~ r~^-^, faster 
than the neutrino density. Moreover, also the neutrino 
densities for tpb < 0.2 s are roughly a factor ^ 3 smaller 



than in the case of the iron-core supernovae. This reflects 
the practical absence of an extended accretion phase for 
this lovif-mass star. In this case the explosion succeeds 
very shortly after the core-bounce. Therefore, the shock- 
front is already beyond the radial range interesting for 
the flavor conversions. 

Fig. [T^ shows the time evolution of i?min in the range 
T G [''sync; ?'cnd] relevant for collective conversions. We 
realize that the behavior of i?min is different with respect 
to the previously considered iron-core supernovae. Since 
in this case the matter does not present an abrupt dis- 
continuity, but it monotonically decreases, the position of 
Rmin will always coincide with rend- Moreover, since the 
electron density is never much larger than the neutrino 
density in the conversion region, it results in i?niin < 3. 

The peculiar time evolution of i?min in this case will 
depend on the small differences between the matter and 
the neutrino term. At tph < 0.13 s, the electron density is 
smaller than the neutrino density in the range of collec- 
tive conversions, giving i?niin<l- Then, the ratio i^min 
grows up to ^ 2 — 3 at tpb — 0.2 s since the electron 
density in the conversion region is larger than in the pre- 
vious case. As a result, the time-evolution of the neutrino 
survival probabilities presents a change between a regime 
dominated by the v-v effects at very early times and one 
of matter-suppressed oscillations at later times. 

In Fig.[T3]we explicitly show this behavior, represent- 
ing the radial evolution of the Vg survival probability Pee 
for the same post-bounce times snapshots of Fig.[TTl with 
(continuous curves) and without (dashed curves) matter. 
We find that at ipb = 0-08 s where i?<l, the matter 
suppression is relatively small and the neutrino-neutrino 
interactions produce a quite almost complete swap be- 
tween Ve and Vx spectra {Pee = 0.15 at the end of the 
evolution). Conversely, for the later times considered, 
the flavor conversions are strongly suppressed with a fi- 
nal Pee — 0.7 — 0.9, since i?> 1- We want to remark that 
the time evolution of the neutrino oscillation probability 
during the accretion phase is significantly different with 
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respect to the case of the iron-core supernovae we stud- 
ied. This will allow the distinction of a 0-Ne-Mg-core 
SN from a iron-core SN, in the case of the detection of a 
future galactic event. 
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FIG. 13: 8.8 Mq progenitor mass. Radial evolution of the 
survival probability Pee for electron antineutrinos at different 
post-bounce times for the multi-angle simulations in presence 
of matter effects (continuous curve) and for the case with 
Tie = (dash-dotted curve). 



EARTH MATTER EFFECT DURING THE 
ACCRETION PHASE 



It is plausible that a high-statistics neutrino detection 
from a future galactic SN [s^, would faithfully moni- 
tor the abrupt time variations imprinted on the neutrino 
oscillation probabilities by the transitions between stages 
of complete and partial matter suppression during the 
accretion phase. We plan to perform a dedicated inves- 
tigation of the matter suppression during the accretion 
phase on the observable SN neutrino burst in a future 
work. 

Here we discuss another consequence of the SN dense 
matter effects, namely the change of the interpretation 
of the Earth matter effect on the supernova neutrino sig- 
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nal. It is widely known that if SN neutrinos arrive at 
the detector from "below" , the Earth crossing will induce 
an energy-dependent modulation in the neutrino survival 
probability 5y] . The appearance of the Earth effect de- 
pends on the neutrino fluxes and on the mixing scenario. 
The accretion phase is particularly promising to detect 
Earth crossing signatures because the absolute SN v flux 
is large and the flavor-dependent flux differences are also 
large. 

Before the inclusion of the collective effects in the SN 
neutrino flavor transitions, the detection of Earth matter 
modulations has been proposed as a tool to distinguish 
the neutrino mass hierarchy at "lar ge" values of the mix- 
ing angle ^13 (i.e. sin^ ^13 > 10"'') [60|. Then, when col- 
lective oscillations were taken into account, it was specu- 
lated that due to the self-induced flavor transformations, 
the presence or absence of Earth matter effects will break 
the degeneracy between normal and inverted mass hi- 
erarchy for "small" values of the mixing angle ^13 (i.e. 
sin^ 013 < 10~^) while no hierarchy discrimination is pos- 
sible for larger values of ^13 [6l| . Now, due to the matter 
suppression effect this conclusion should be revisited for 
the accretion phase. 

To determine the observable SN neutrino fluxes at 
Earth, we refer to [s^l. We consider the modified flavor 
basis, {ve,Vx,Vy), which is defined such that [vejVxjVy) = 

-R23^(^23)(^'ej '^/ij '^r); whcrc i?23 IS the 2-3 rotation ma- 
trix. This will allow for collective oscillations between e 
and y states. For definiteness, here we concentrate on 
the Vf, spectra observable through inverse beta decay re- 
actions, Ve + p — > n + e"*" , at water/ice Cherenkov or 
scintillation detectors [53, [62] . 

In normal hierarchy, the Fe flux at Earth, F-P , (without 
collective effect and no MSW effect associated to ^13) for 
any value of the mixing angle ^13 is given by (60j 

= cos^ 0i2Fp^ + sin^ 0i2Fp^ , (16) 

where 6*12 is the 1-2 mixing angle, with sin^ 612 ~ 0.3 [sl. 

In the inverted mass hierarchy, MSW matter effects in 
the SN envelope are characterized in terms of the level- 
crossing probability Ph of antineutrinos, which is in gen- 
eral a function of the neutrino energy and ^13 [36j . In 
the following, we consider two extreme limits, namely 
Ph when sin^ 6*13 > lO'^ ("large"), and Ph 1 
when sin^ 6*13 < 10~^ ("small"). In this situation and for 
complete matter suppression of collective oscillations for 
large 6*13 , F^^ is given by 

Fi = F,^ , (17) 

while small 613 results in the same case as for normal 
hierarchy [Eq. (fT6)) ]. 

In the case of matter-induced docoherence between i7g 
and Vy in inverted hierarchy, leading to a complete mix- 
ture of T^e and Vy, the detectable flux for any value of 
the mixing angle 6*13 would be given by 

F£ = cos^ 9,2 ^'^^/'" + sin^ e^2F,^ . (18) 
z 



Earth effect can be taken into account by just mapping 
cos^ 012 — P(l7i — Ve) and sin^ 0i2 — )■ 1 — P{vi — Ve). 
where P{vi ~->Ve) is the probability that a state entering 
the Earth as mass eigenstate Vi is detected as Ve at the 
detector [59| . 

Considering for definiteness the case of iron-core SNe, 
the presence or absence of Earth matter effects at early 
times (ipb < 0.2 s) will allow to distinguish the neutrino 
mass hierarchy at large value of the mixing angle ^13 
[see Eqs. p^ - ((T7)) ]. Moreover, the Earth effects may 
become observable via the interesting time pattern in- 
troduced above, in the case of inverted mass hierarchy 
at large ^13: (1) no effects at early times [Eq. ([T7|]. (2) 
the appearance at intermediate times when decoherence 
occurs [Eq. ((T5)) ]. and (3) disappearance again when the 
accretion rate decreases significantly. 



VI. DISCUSSION AND CONCLUSIONS 

Different simulations of core-collapse SNe [sol - li^ . 
show independently that the matter density in the stellar 
interior is large during the post-bounce accretion phase 
before the onset of an explosion. This implies that self- 
induced neutrino flavor transformations are affected by 
the high matter density, through trajectory-dependent 
phenomena [37|. 

In order to characterize the SN ly flavor evolution quan- 
titatively, we performed a dedicated study, where we 
took as benchmark for the SN neutrino emissivity and 
the matter proflles the recent long-term core-collapse SN 
simulations from Ref. We considered three different 
cases for the supernova progenitor mass, the low mass 
8.8 Mq 0-Ne-Mg-core and the two more massive iron- 
core progenitors of 10.8 and 18 Mq. In all three cases the 
electron density Ue is never negligible in comparison to 
the neutrino density n„ during the accretion phase. As a 
consequence, the trajectory-dependent matter effects al- 
ways influence the development of the self-induced trans- 
formations. We realized that during the accretion phase 
both the condition of matter suppression (ue 3> n.^) and 
matter-induced decoherence {ui, ~ rig) are realized. This 
implies that interesting time-dependent features in the 
neutrino signal can occur at early times. Moreover, the 
different patterns of complete/partial matter suppression 
would allow, at least in principle, to distinguish from 
the flavor evolution between iron-core and 0-Ne-Mg core 
SNe. 

In early schematic investigations, the accretion phase 
seemed a particularly clear setup to probe the devel- 
opment of the collective neutrino transformations, be- 
cause the expected flux hierarchy is large and robust. In 
such a situation the pattern of the self-induced spectral 
swaps/splits were thought to be unambiguous |32l |. In 
this scenario, the presence or absence of Earth matter 
effects on the SN neutrino burst were proposed to break 
the degeneracy between normal and inverted mass hi- 
erarchy for "small" values of the mixing angle ^13 [6l| . 
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Our findings seem to change this picture once more. In 
particular, the dense matter suppression of cohective os- 
cillations in supernovae implies that Earth matter effects 
may allow to distinguish the neutrino mass hierarchy in 
the likely case of large 9i3 _6J|] , as originally expected in 
the analysis performed without the inclusion of the col- 
lective phenomena. This intriguing possibility makes su- 
pernova neutrino detection a competitive way to get this 
yet unknown property of the neutrino mass spectrum in 
addition to terrestrial experiments [s^l . 

Furthermore, in the past [H, Ull it has been speculated 
that swaps of neutrino fluxes between electron and non- 
electron flavors in the deepest SN regions may increase 
the neutrino energy deposition and hence the neutrino 
heating in the gain region behind the standing bounce 
shock. It was argued that it could possibly help to re- 
vive the shock and to trigger neutrino-driven explosions. 
We have shown that the presence of the matter suppres- 
sion of flavor conversions at high densities, behind the 
shock front (see, e.g.. Fig. 6), implies that these cannot 
play any significant role in changing the neutrino energy 
deposition rate behind the stalled shock-wave. As an im- 
portant consequence, the problem of he v flavor mixing 
in SNe can be decoupled from the v transport and im- 
pact on the matter heating/cooling. This result is con- 
flrmed by the analysis recently performed in [65| using 
two-dimensional supernova models. In that work, even 
if multi-angle matter effects are not included in the sim- 
ulation of the flavor oscillations, a negligible impact of 
the self-induced conversions is found on the energy depo- 
sition behind the shock-wave. Indeed, we also find that 
neglecting matter effects, the possible range where self- 
induced oscillations would have significantly developed is 
always after the shock- wave position (rgh^'^sync) at the 
times relevant for the shock revival (see our Fig. 6-8). 

Our results have been obtained considering a spheri- 
cally symmetric neutrino emission. All the previous anal- 
ysis in the field have relied on this assumption to make 
the flavor evolution equations numerically tractable. It 
remains to be investigated if the removal of a perfect 
spherical symmetry can provide a different behavior in 
the flavor evolution [66| . Moreover, in multi-dimensional 
SN models density fluctuations are expected behind the 
standing bounce shock, due to the presence of convection 
and hydro instabilities. These can range at most between 
10% to a factor 2-3 (see, e.g., [iO,!!^)- Therefore, even in 



the case a fraction of neutrinos beam crosses underdenses 
regions, the matter suppression of the collective oscilla- 
tions will still remain relevant. This claim is supported 
by a recent analysis of the matter suppression, performed 
with two-dimensional SN simulations 65] . 

Self-induced effects in the supernova neutrinos still re- 
main crucial during the later cooling phase, when the 
matter density decreases continuously {tph > 1 s) and be- 
coming sub-dominant with respect to the neutrino den- 
sity. In this situation and for the flux ordering provided 
by the simulations from Ref. Q, we conflrm that self- 
induced oscillations in the cooling phase produce multi- 
ple splits and swaps, like the one described in [26l - [3^ . 
Moreover, for low-mass 0-Ne-Mg-core SN a peculiar in- 
terplay of MSW and self-induced effects would produce 
interesting signatures during the i>e prompt neutroniza- 
tion burst |68l470l |. 

In conclusion, the time evolution of the SN neutrino 
signal provides a large amount of information on the fla- 
vor transformations of neutrinos in the deepest stellar 
regions. The early post-bounce deleptonization, accre- 
tion and cooling phases represent three different "exper- 
iments" that allow us to probe different behaviors of the 
dense SN neutrino gas. In order to learn the most from 
a future neutrino observation from a galactic event, im- 
proved theoretical and experimental studies are necessary 
to decipher the supernova neutrino enigma. 
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